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ABSTRACT 

(N 
O 

^ I An analytic model of the evolution of a rotating black hole (BH) is proposed by consid- 

p • ering the coexistence of disk accretion with the Blandford-Znajek process. The evolutionary 
Q . characteristics of the BH are described in terms of three parameters: the BH spin a^,, the 
^ [ ratio k of the angular velocity of the magnetic field lines to that of the BH horizon and the 
' parameter A indicating the position of the inner edge of the disk. It is shown that the ratio k 
being a little greater than 0.5 affects the evolutionary characteristics of the BH significantly, 
and the BH spin increases rather than decreases in its evolutionary process provided that 
^ I the initial value of the BH spin is located in an appropriate value range determined by ratio 
k. Our calculations show that the system of a BH accretion disk with k = 0.6 might provide 
' a much higher output energy in a shorter timescale for gamma-ray bursts than the same 
system with k = 0.5. 

Subject headings: Black hole — accretion disk — magnetic fields 
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^ ! system with k = 0.5. 
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^ : 1. INTRODUCTION 

As is well known, the Blandford-Znajek(BZ) process was proposed originally as a possible energy 
mechanism of quasars and active galactic nuclei(AGNs;Blandford & Znajek 1977; Macdonald & Thorne 
^ I 1982, hereafter MT82; Rees 1984). Recently, Lee, Wijers, & Brown (2000) proposed that the BZ process 
■ can be used as a central engine for powering gamma-ray bursts (GRBs), where the rotating energy of 
a stellar black hole (BH) with a magnetic field of lO^^G is extracted along the magnetic field lines 
supported by a transient magnetized accretion disk. Very recently, Lee and Kim (2000, 2002) proposed 
a model of the evolution of a rotating BH at the center of GRBs by considering the effects of the BZ 
process with the transient disk (hereafter the LK model), and some constraints to the parameters of 
the BH accretion disk are given. 

In this paper an analytic model of BH evolution is proposed based on the LK model and is improved 
in two respects: (1) A parameter A is introduced for the position of the inner edge of the disk/torus, 
which is located between the innermost bound orbit and the last stable orbit. (2) Another parameter k 
is used to indicate the ratio of the angular velocity of the magnetic field lines to that of the BH horizon. 
It is shown that the ratio k being a little greater than 0.5 affects the evolutionary characteristics of the 
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BH significantly, and the BH spin increases rather than decreases in its evolutionary process provided 
that the initial value of the BH spin is located in an appropriate value range determined by the ratio 
k. Our calculations show that the system of a BH accretion disk with k — 0.6 might provide a much 
higher output energy in a shorter timescale for GRBs than the same system with k — 0.5. 

This paper is organized as follows: In §2 wc derive the basic equations of BH evolution and the 
corresponding characteristic functions by introducing two parameters A and k. In §§3 and 4 we discuss 
the effects of the parameters k and A on BH evolution and GRBs by using the parameter space and the 
corresponding characteristic functions, respectively. Finally, in §5 we discuss some problems concerning 
our model. Geometrized units {G — c— 1) are used in this paper. 



2. BASIC EQUATIONS OF BH EVOLUTION AND CHARACTERISTIC 

FUNCTIONS 

Now we are going to describe our model for the evolution of a rotating BH surrounded by a 
magnetized accretion disk based on the LK model. First, a parameter A is introduced to indicate the 
position of the inner edge of the disk in the following equation: 

Xin^Xmb + HXms-Xmb): < A < 1, (l) 

where Xin is the dimensionless radial parameter corresponding to the inner-edge radius rj„ of the disk, 
which is located between rmb and r^s (Abramowicz, Jaroszynski, & Sikora 1978; Abramowicz & Lasota 
1980). The quantities 

Xin: Xmhi ^-ud Xms ^-re dimensionless radial parameters: 



^'""V^' ^""^^y^' ^""'^yw 

It is found from equation (1) that the inner edge of a thick disk varies continuously from the innermost 
bound orbit to the last stable orbit as A changes from zero to unity, and the disk becomes a thin disk 
as A attains unity. 

Usually the evolutionary state of a Kerr BH can be described by two parameters: BH mass M and 
spin a*, and the latter is related to the BH mass M and angular momentum J by 

a* = 0<a^<l. (3) 

Both Xmb and Xms depend only on a* (Novikov & Thorne 1973; Abramowicz et al. 1978): 

Xmb = 1 + VI - a*, (4) 

Xts - 6xL + 8a^Xms - 3al = 0. (5) 

The specific energy Ein and the specific angular momentum Lin corresponding to the inner-edge radius 
Tin of the disk/torus are expressed by (Novikov & Thorne 1973) 
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In the LK model, the effects of the BZ process on the evolution of a rotating BH are taken into account 
in the environment of a transient magnetized accretion disk. Due to the lack of knowledge of the remote 
astrophysical load in the BZ process, it is doubtful that the matching condition k = 0.5 can be satisfied 
exactly. Therefore, we introduce a parameter k to describe the rotation of the magnetic field lines 
relative to the spinning BH, 

k = ^, 0<k<l, (8) 
"ft 

where Qf is the angular velocity of the magnetic field lines and is that of the horizon and is related 
to the horizon radius by 

nn^^, r, = M{l + q). (9) 

h 

By using a modified equivalent circuit for the BZ process based on MT82, we derive the expression for 
the rate of extracting energy from the rotating BH in the BZ process (hereafter the BZ power; Wang, 
Xiao, & Lei 2002) as follows: 

Pbz = 4A;(1 - A;)P^f (10) 
po,umai ^ BlM'lQ-'arctanQ - i^], (11) 



g=^^, (12) 

where p^^*"*"' jg ^j^e optimal BZ power corresponding to the impedance matching with k = 0.5 and Bh 
is the magnetic field on the horizon. Both Q and q are the parameters depending only on the BH spin. 
Taking k = 0.5 in equation (10), wc have the same expression for the optimal BZ power as given by 
Lcc ct al. (2000), who pointed out that the BZ power had been underestimated by a factor of 10 in 
previous works. 

Based on the conservation of energy and angular momentum, the basic evolutionary equations of a 
rotating BH with disk accretion by the BZ process can be written as (Park & Vishniac 1988; Moderski 
& Sikora 1996; Moderski, Sikora, & Lasota 1997; Wang, Lu, & Yang 1998; Lee & Kim 2000) 

^ = Eir^Md - Pbz: (13) 

Incorporating equations (13) and (14), we obtain the evolutionary equation of the BH spin: 

^ = M-\UnMa - ^) - 2M-'a4E,nMd - Pbz). (15) 
at K\lh 

The accretion rate of the transient magnetized disk is related to Bh by (Lee & Kim 2000) 
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where w is the Kerr metric parameter at the inner edge. According to the LK model, depends on 
the mass loss from the disk (Lee & Kim 2000): 

Bl = Bl{^)D{t), (17) 

\l Mddt , , 

where Md(0) and -B?i(0) are the initial disk mass and magnetic field, respectively. 

Substituting the above expression for Ein^ Lin, Xin, ^h, Pbz and into equations (13) and (15), 
we obtain the following evolutionary equations: 

^^f{a^,k,X)M,, (19) 

^ = M-'g{a.,k,X)M,. (20) 

From equations (19) and (20), we find that dM/dt and da^^/dt have the same sign as f{a^,k,X) and 
g{a*,k,X), and hereafter the two functions are referred to as the characteristic function of BH mass 
(CFBHM) and that of BH spin (CFBHS), respectively. 



3. EFFECTS OF THE PARAMETER k ON BH EVOLUTION AND GRBs 

The characteristics of BH evolution can be well described by using the CFBHM and CFBHS, and 
the curves of these two functions varying with a* for a thin disk (A = 1) with different values of k 
are shown in Figure 1. It is shown in Figure la that /(a*,A;, 1) is negative when a* > 0.9981,0.9984, 
and 0.9999 for k = 0.5, 0.6, and 0.8, respectively. From Figure lb we find that g{a^, k, 1) is always 
negative for k — 0.5, while it is positive for k — 0.6 with 0.3114 < a* < 0.8067, and for k — 0.8 with 
0.1205 < a^, < 0.9598. These results imply that the BH mass might decrease as the BH spin approaches 
unity very closely, and the BH could be spun up in the duration of powering the GRB, provided that 
the parameter A; is a little greater than 0.5. 

The effect of the variation of k on the evolution of the BH can be discussed more visually in the 
parameter space consisting of a* and k with the curves represented by /(a*, k,l) = and g{a^, k, 1) = 0, 
as shown in Figure 2. 

It is found that the parameter space is divided into three regions by two thick solid curves, which 
are represented by 9^(0*, k,l) — and /(a*, k, 1) = 0. The former is the common boundary of regions I 
and II (curve abc), as shown in Figure 2a, and the latter is the common boundary of regions I and HI, 
as shown in Figure 2b. Regions I and HI are each further divided into two sub-regions, lA and IB, and 
IIIA and IIIB, by the horizontal thin line bd, as shown in Figures 2a and b, respectively. The segment 
bd is tangent to the curve g{a^,,k,l) = at the bottom point b(0.5610, 0.5118). In the parameter 
space, each filled circle with an arrow is referred to as a representative point (RP), which represents 
one evolutionary state of BH. We can use RP displacement in the parameter space to describe the 
evolution of the BH. According to the positions of RPs, we have five different BH evolutionary states 
corresponding to five subregions, as shown in Table 1. 
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From Table 1 we find that the five evolutionary states of a rotating BH can be easily determined 
by RP position in the parameter space, and the details are given as follows: 

(1) The RPs in regions lA and IB represent the evolutionary states of the BH with increasing mass 
and decreasing spin, and the BH with the RP in region lA will never stop rotating, while the BH with 
the RP in region IB might evolve to a Schwarzschild BH with zero spin in long enough time. The 
evolution state given in the LK model is just described by the RP in region IB. 

(2) The RPs in region II represent the evolutionary states of the BH towards the equilibrium spin 
with increasing mass and spin increasing toward equilibrium, which makes the BH a more powerful 
central engine for GRBs. 

(3) The RPs in regions lllA and IlIB represent the evolutionary states of the BH with decreasing 
mass and decreasing spin, which implies that the efficiency of transferring the accreted mass into GRB 
energy might be greater than unity in these two evolutionary states. 

(4) As shown in Figure 2 and Table 1, there are two possible evolutionary terminals for the BH, if 
the evolution time is long enough. The terminal with equilibrium spin a^* corresponds to the RPs in 
regions lA, II, and IIIA, while that with zero spin corresponds to the RPs in regions IB and IIIB. 

Prom the above discussion, we find that evolutionary characteristics of a rotating BH in the BZ 
process are very sensitive to the parameter k, especially as k varies across the vertical ordinate of the 
bottom point ~ 0.5118 in Figure 2a. 

The value range of a* in region II is expressed by a*"'"" < a* < a^^, where a*"™ and al'^ arc, 
respectively, the left and right intersections of the curve abc with a horizontal line above the segment 
bd. Inspecting Figure 2a, we find that the more the horizontal line is above the segment bd or the 
greater k is than k^°*, the greater the value range of a* is in region II. The left part of curve abc is the 
set of a*"*"", while the right part of curve abc is the set of a^*. Obviously, a*"^" is unstable, while a^* is 
stable for the evolutionary states. 

In order to discuss the effects of the parameter k on the evolution of a rotating BH and the related 
characteristics of GRBs, we compare the results of the BH accretion disk system by taking k = 0.6 and 
k = 0.5, respectively. Our main results are shown in Figures 3-8, in which the solid lines are used for the 
results with k — 0.6, and the dashed hues for those with k — 0.5. In our calculations, the initial values 
of the concerned quantities are taken as those given in the LK model: M/i(0) = 7Mq, Md{0) — SM©, 
and Bd{0) = lO^^G. 

(1) In Figure 3 the BH spin a* increases with the evolution time if its initial value is taken to be 
0.3114 < a*(0) < 0.8067 for k = 0.6, while it always decreases with the evolution time no matter what 
initial value it takes for k = 0.5. As shown in Figure lb, 0.3114 and 0.8067 are exactly the values of 
Qturn Q^^^ a^* for = 0.6, respectively. 

(2) In Figure 4 can never be greater than a*"*^'" = 0.3114 no matter how much the initial disk 
mass is for k — 0.6, while it increases monotonically with the initial disk mass for k — 0.5. Since al is 
defined as the critical BH spin to ensure that the BH will stop spinning after all the initial disk mass 
accretes onto the BH in the LK model, our result implies that a BH with its initial spin greater than 
Qturn never be spun down to a Schwarzschild BH. 

(3) In Figure 5 the output energy E from a system with k = 0.6 {solid lines) is always greater 
than that from a system with k = 0.5 {dotted lines). In the former case, E will monotonically increase 
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with the initial disk mass Md(0) without an upper hmit, provided that the initial BH spin is greater 
than a*"^", while in the latter case, it always keeps constant for Ma{0) greater than one critical value. 
Another difference between these two systems lies in the derivative of E with respect to Md{0). As 
shown in Table 2 the first derivative dE/dMd{0) always increases with Mdi^O) (in units of solar mass) 
provided that the system has the parameter k = 0.6 with 0.3114 < a*(0) < 0.8067, while it always 
decreases with Md(0) either for k — 0.6 with < a*(0) < 0.3114 or for k — 0.5 with any value of a*(0). 
These results can be easily explained by the RPs in the parameter space (Figure 2a) and the relation 
of the BZ power to the BH spin. 

(4) In Figure 6 the output energy from the system with k = 0.6 is always greater than that with 
k = 0.5 for any initial values of the BH spin. Defining the relative difference of output energy between 
the above two systems as S = {E^^^ — E^^)/E^^, we have the curve of 5 versus a*(0) as shown in 
Figure 7, and the maximum is 5jnax ~ 2.44 at a*(0) 0.294. Thus, the output energy for GRBs can be 
remarkably augmented by a little increase above k^°*. 

(5) In Figure 8 the evolutionary timescales tevi of GRBs versus a*(0) are obtained by incorporating 
equations (16) — (18) and the initial values of the systems with = 0.5 and = 0.6, where the cutoff 

of tevi is taken as Tgo defined in the LK model. A peak value of t^vi is found for each system, which is 
located at a*(0) ~ 0.298 for k = 0.6 and at a^(0) ~ 0.402 for k = 0.5, respectively. It is shown that the 
tevi of the system with k = 0.6 is always less than that with k = 0.5 for a*(0) > 0.366. These results 
might be helpful to explain GRBs with the most energetic power. 



4. EFFECTS OF THE PARAMETER A ON BH EVOLUTION AND GRBs 

Taking the parameter A into account, we can discuss the evolutionary characteristics of a BH 
surrounded by a thick disk/torus, and the parameter spaces with boundary curves g{a^, k, \) = and 
f{a^,k,\) = corresponding to different values of A are shown in Figure 9. It is found that the 
above five evolutionary states of a BH surrounded by a thin disk (A = 1) remain valid in the case of 
a thick disk/torus with its inner edge positioned by the parameter A (0 < A < 1). It is interesting 
to note that both the bottom point of the curve g{a^, k,X) — and the leftmost point of the curve 
f{a^,k,\) = vary nonmonotonically with the parameter A, and we have the curves of k^°^ and a^*^-^* 
versus A as shown in Figures 10 and 11, respectively. The quantity A:^°* attains its maximum of 0.598 
at A ~ 0.137, and a'/-^* attains its minimum of 0.995 at A ~ 0.261. The output energy E^ut from the 
systems and the evolutionary timescales tevi corresponding to the different values of A with k — 0.6 are 
compared in Figuresl2 and 13, respectively. Finally, the relative differences of output energy defined 
hj 6 = (-Eo„( — E^'^^)/E^'^^ corresponding to the different values of A are shown in Figure 14. The main 
effects of the parameter A on BH evolution and GRBs arc summarized as follows. 

(1) The BH will be spun up in the case of a thick disk (0 < A < 1) as well as in the case of a thin 
disk (A = 1), provided that we have k > k^°^ with a*"*""" < a*(0) < a^'^, where a*"*"" and a^^ have the 
same meaning as given in §3. 

(2) The efficiency of transferring the accreted mass into GRB energy might be greater than unity 
provided that the RP of BH evolutionary state is located in the corresponding region HI as shown in 
Figure 9b. 

(3) From Figures 12 and 13 we find that the output energy from a system with a thick disk/torus 
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might be greater, and the concerned evolutionary timescales might be longer, than the corresponding 
values of a system with a thin disk if a*(0) is greater than some value. Therefore, our model with 
different values of A might adapt to the GRBs with different output energy and different timescales. 

(4) Prom Figure 14 we find that the output energy from a system with k — 0.6 is always greater 
than that from a system with k — 0.5, no matter what value a*(0) has. This is true for a system with 
a thick disk/ torus as well as for a system with a thin disk. 

In summary, compared to the effects of k the characteristics of BH evolution and GRBs are not 
remarkably affected by variation of the parameter A. 



5. DISCUSSION 

In this paper the evolutionary characteristics of a rotating BH surrounded by a transient magnetized 
disk are discussed by considering the coexistence of disk accretion with the BZ process, and the two 
parameters A and k are introduced in our model to modify the LK model. A main consequence of our 
model is that the BH will be spun up rather than spun down provided that k is greater than a critical 
value with " < O'*(0) < a^''. Our calculations show that more energy from a BH system with k 
above A;^"* might be extracted in less time compared to a BH system with k below k^°*. 

However there are still some problems related to our model. First, the different values of the ratio 
k = Qf/Qh play a very important role in our model. Unfortunately, the determination of this ratio 
has remained one of the main unresolved problems in the BZ process. It was argued in a speculative 
way that the ratio might be regulated to about 0.5, if the charged particles conspired with the BH to 
have just the right amount of inertia for the impedance matching (MT82). However, it is difficult to 
understand how the load can conspire with the BH to have the same resistance and satisfy the matching 
condition, since the load is so far from the BH that it cannot be causally connected (Punsly & Coroniti 
1990). 

It was pointed out in the theory of BH magnetospheres that flf can be regarded as a function 
of magnetic flux ^, i.e., since it is constant on the magnetic surface due to isorotation and 

axial symmetry (MT82). Very recently, Beskin and Kuznctsova (2000) discussed the stream equation 
describing magnetic surfaces v[f(r, ^), and found that Qf can be either greater or less than for 
the different possibilities of the flow, with the electric current fixed by the pair creation region in the 
BH magnetosphere. However, we still have a long way to go for the determination of Q/ in a BH 
magnetosphere. 

Second, we notice that the relation in equation(16) between and is the key to affecting the 
evolutionary characteristics. To illustrate this, we replace equation (16) with another relation between 
Md and B^, which was proposed based on the balance between the pressure of the magnetic field on 
the horizon and the ram pressure of the innermost parts of the accretion flow (Moderskin et al. 1997), 
i.e., 

^ = Pram ^ ^ (21) 

From equation (21) we assume the relation to be 

M, = V. (22) 
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Hereafter equation (22) is referred to as the MSL relation. Substituting the MSL relation into equation 
(13) and (15), we obtain the following evolutionary equations: 

^-fMsA(^*,k,X)Ma, (23) 

^^M-'g^,,ia.,k,X)m, (24) 

where f msA'^*-i^-i^) ^^"^ 9msA'^*-i^-i^) CFBHM and CFBHS corresponding to the MSL relation, 
respectively. In the same way, we obtain the parameter space consisting of a* and k as shown in Figure 
15. The parameter space is very different from that depicted in Figure 2, and it is only divided into 
two regions by the boundary curve 5'j,,^5^(a*, k, A) = 0, corresponding to only two states of BH evolution. 
The characteristics of BH evolution related to the MSL relation are shown in Tabel 3. 

In this paper an analytic model of BH evolution is proposed for GRBs. Although this model is 
rather particular and speculative , the potential application is still attractive: a little variation of the 
angular velocity of the magnetic field lines might affect the output energy and evolutionary timescales 
of a BH system remarkably. Wc also hope to improve and extend this model to AGNs and stellar-mass 
BHs in future work, since there is more evidence to strongly suggest the existence of fast-spinning BHs 
and possibly energy extraction via Blandford-Znajek-like processes(Krolik 2001). 

The anonymous referees are thanked for their suggestions on the improvement of our manuscript. 
V. S. Beskin is thanked for his encouragement and the discussion of the determination of f2j in a BH 
magnetosphere. This work is supported by the National Natural Science Foundation of China under 
grant No. 10173004. 
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Table 1. BH EVOLUTIONARY STATES CORRESPONDING TO FIVE DIFFERENT 

SUBREGIONS 



Region CFBHM CFBHS RP Displacement Evolutionary States Terminal 



lA > < towards the left . . equilibrium spin 

mcreasmg mass, 



Tr> n n X J XI 1 decreasing spin 

IB > < towards the leit zero spm 



TT r. n X 1 XI • 1 J. mcreasmg mass, . 

11 > U > U towards the right . . . equilibrmm spm 

mcreasmg spm 



IIIA < < towards the left , . equihbrium spin 

decreasmg mass. 



TTTT^ n n X J XT T rx decreasing spin 

lllB < U < U towards the left zero spm 
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Table 2. dE/dMd{0) VARYING WITH Md{0) AND a,(0) FOR THE TWO SYSTEMS 



a*(0) k 



0.05 0.10 0.15 0.20 0.25 



0.2 0.5 0.1868 0.1849 0.1826 0.1796 0.1749 

0.6 0.1827 0.1818 0.1808 0.1797 0.1784 

0.3 0.5 0.1996 0.1989 0.1981 0.1974 0.1965 

0.6 0.19495 0.19489 0.19484 0.19478 0.19473 

0.4 0.5 0.2135 0.2133 0.2130 0.2127 0.2124 

0.6 0.2085 0.2088 0.2091 0.2094 0.2097 

0.5 0.5 0.2299 0.2298 0.2297 0.2296 0.2295 

0.6 0.2245 0.2250 0.2254 0.2259 0.2263 

0.6 0.5 0.2502 0.2501 0.2501 0.2500 0.2499 

0.6 0.2445 0.2450 0.2455 0.2460 0.2465 

0.8 0.5 0.3186 0.3175 0.3165 0.3155 0.3145 

0.6 0.31187 0.31193 0.31198 0.31203 0.31208 
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Table 3. BH EVOLUTIONARY STATES CORRESPONDING TO THE MSL RELATION 



Region 


CFBHM 


CFBHS 


RP Displacement 


Evolutionary States 


Terminal 


I 


> 


< 


towards the left 


increasing mass, 
decreasing spin 


equilibrium spin 


II 


> 


> 


towards the right 


increasing mass, 
increasing spin 


equilibrium spin 
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Fig. 1. — (a) Curves of /(a*, /c, A) vs. a*, which is referred to as CFBHM to indicate the evolutionary 
characteristics of the BH mass, (b) Curves of g{a^,k,\) versus a*, which is referred to as CFBHS to 
indicate the evolutionary characteristics of the BH spin. The parameters for the model are A = 1, 
k — 0.5 {solid line), k — 0.6 {dashed line) and k — 0.8 {dotted line). 
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Fig. 2. — Parameter space for BH evolution: the boundary curve of regions I and II is represented 
by g{a^,, k, A) = 0, and the boundary curve of regions I and III is represented by /(a*, k, A) = 0. BH 
evolutionary state is represented by the filled circles with arrows (RPs). It is noted particularly that the 
RP in region II moves towards the right, which means the BH spin is increasing rather than decreasing 
in its evolutionary process. The parameters for the model are A = 1 and < k < 1. (a) < a* < 1; (b) 
0.997 < a* < 1 . 



- 15 - 

0.9 
0.8067 

0.7 

0.6 

0.5 

0.4 
0.3114 

0.2 
0.1 
0.0 

1000 2000 3000 4000 

t[s] 

Fig. 3. — Evolution of the BH spin with time in the units of seconds with sohd hues for k = 0.6 and 
dashed hues for k = 0.5. It is noted that the BH spin is increasing with time as shown by the sohd hues 
with 0.3114 < a,(0) < 0.8067. The parameters for the model are M/,(0) = 7Mq, Md(0) = 3Mq and 




0.6 F 



0.5 - 



0.4 - 



* 0.3 



0.2 - 



0.1 - 



0.0 y 








10 20 30 40 



50 



Fig. 4. — Critical BH spin to ensure that the BH will stop spinning after all the initial disk mass 
accretes onto the BH. It is noted that the solid line {k = 0.6) is always below al = 0.3114 for any values 
of Mrf(O). The Parameters for the model are the same as for Fig. 3. 
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Fig. 5. — Output energy vs. the initial disk mass with various initial BH spin a*(0). It is noted that 
the solid line {k — 0.6) is always higher than the dashed line(A; = 0.5), and the former has no upper 
limit for a*(0) > 0.3114. The parameters for the model are the same as for Fig. 3 with a*(0) taken as 
(a) 0.2, (b) 0.3, (c) 0.4, (d) 0.5, (e) 0.6, (f) 0.8. 
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Fig. 6. — Upper limits of the energy out of the system vs. the initial BH spin. It is noted that the solid 
line {k = 0.6) is always higher than the dashed line(/c = 0.5). The parameters for the model are the 
same as for Fig.3. 
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Fig. 7. — Relative difference of output energy vs. the initial BH spin with its 2.44 at 

a*(0) ~ 0.294. The parameters for the model are the same as for Fig.3. 
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Fig. 8. — Evolutionary timescale of a GRB powered by two BH systems with k = 0.6 {solid line) 
and k = 0.5 {dashed line). It is noted that the former is less than the latter for a*(0) > 0.366. The 
parameters for the model are the same as for Fig. 3 with the cutoff of tg^i taken as Tgo in the LK model. 
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Fig. 9. — Parameter space for BH evolution corresponding to different values of A. It is noted that the 
boundary curves vary nonmonotonically with A . The parameters for the model are k = 0.6, A = 1 
{solid line), A = 0.5 {dashed line), and A = {dotted line) with (a) < a* < 1, (b) 0.9956 < a* < 1. 
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Fig. 11. — Horizontal ordinate ajf-^* of the leftmost point of CFBHM varying with A for < A < 1. 
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Fig. 12. — Upper limits of the energy out of the system vs. the initial BH spin. It is noted that the 
output energy from a system of a slow-spinning BH with a thin disk might be greater than that from a 
system with a thick disk, while the output energy from a system of a fast-spinning BH with a thin disk 
might be less than that from a system with a thick disk. The parameters for the model are k — 0.6, 
A = 1 {solid line), A = 0.5 {dashed line), and A = {dotted line). 
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Fig. 13. — Evolution time of a GRB vs. the initial BH spin. It is noted that the evolutionary timescale of 
a slow-spinning BH with a thin disk might be longer than that with a thick disk, while the evolutionary 
timescale of a fast-spinning BH with a thin disk might be less than that with a thick disk. The 
parameters for the model are the same as for Fig. 12. 
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Fig. 14. — Relative difference between the output energy from the system with k = 0.6 and that with 
k — 0.5 for different values of A. It is noted that 5 varies nonmonotonically with a*(0) for different peak 
values corresponding to different values of A and a*(0). The parameter for the model are the same as 
for Fig. 12. 




Fig. 15. — Parameter space for BH evolution corresponding to the MSL relation with different values 
of A. The parameters for the model are A = 1 {solid line), A = 0.5 {dashed line), A = {dotted line). 



